Formation of clusters in the ground state of the t — J model on a two leg ladder 



O 

o 

(N 

> 

o 
O 



A. Fledderjohann 1 , A. Langari 2 and K.-H. Mutter 1 
1 Physics Department, University of Wuppertal, 42097 Wuppertal, Germany and 
2 Institute for Advanced Studies in Basic Sciences, Zanjan 45195-159, Iran 
(Dated: February 2, 2008) 

We investigate the ground state properties of the t — J model on a two leg ladder with anisotropic 
couplings (t, a = J/t) along rungs and (t',a' = J' /t') along legs. We have implemented a cluster 
approach based on 4-site plaqettes. In the strong asymmetric cases a /a' <C 1 and a' /a <C 1 the 
ground state energy is well described by plaquette clusters with charges Q — 2,4. The interaction 
between the clusters favours the condensation of plaquettes with maximal charge - a signal for 
phase separation. The dominance of Q = 2 plaquettes explains the emergence of tightly bound 
hole pairs. We have presented the numerical results of exact diagonalization to support our cluster 
approach. 

PACS numbers: 71.10.Fd,71.27.+a,75.10.Jm 
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I. INTRODUCTION: MOTIVATION AND 
QUESTIONS OF INTEREST 

After the discovery of high T c superconductors 1 - al- 
most twenty years ago - models of strongly correlated 
electron systems doped with holes have attracted much 
interest. The t — J model in two dimensions2iiMi£i£iL£i2i 
and on laddersifi*ii has been studied intensively in order 
to understand the behaviour of mobile holes in an anti- 
ferromagnetic background. Although some exact results 
in the special form of interactions exisliiSiiii the ground 
state and low energy excitations on the t — J ladder are 
not known exactly. 

The generic mechanisms which explain the most strik- 
ing features are of special interest, namely: 

i) the opening of a charge transfer gaiA 

ii) the spatial separation of phases with hole rich and 
hole poor domainesi^. 

Signatures of these features can be seen already in the 
ground state energy per site 



develops a plateau with a width 



s(p) = W , 



Qtot 

N • 



(1.1) 



where Eg, N, Qtot are the ground state energy, the num- 
ber of sites and the total charge, respectively. 

i) A discontinuity in the first derivative, i.e. chemical 
potential, 



ds 
dp' 

(i- p -> po - 
M+ P -> Po + 0, 



(1.2) 
(1.3) 



signals the opening of a gap. 
The inverse function 



A 



P+ - P- 



(1.5) 



which is related to the charge transfer gap. This 
is quite analogous to the plateaux in the magneti- 
zation curve M — M(B), which are related to the 
spin gap. 

ii) A linear behaviour of e{p) in some interval, p\ < 

P < P2- 



£l(p) = 



1 



[e(pi)(p 2 -p) + e(p 2 )(p-pi)}, (1.6) 



P2 - Pi 

signals the spatial separation of two phases: the 
first one with charge density pi, the second one 
with p 2 . 

For demonstration, let us consider a two leg lad- 
der, which we divide into two parts with sites 
Nj, charges Qj, Hamiltonians Hj, ground state 
wave functions tp(pj, Nj) and ground state energies 
e(pj)Nj(p) as illustrated in Fig. ^ 



Ni( P ) 



P2 



-N, 



p2 - Pi 

Q 1 =N 1 (p) Pl , 
HpuN^p)), 
e( Pl )Ni(p), 



H 2 (N 2 (p)) 
N 2 ( P ) = 

p2 - Pi 

Q 2 = N 2 (p)p 2 

tp(p 2 ,N 2 (p)) 

£{p 2 )N 2 (p) 



P(p) = Po for p,- < p < p, + , 



(1.4) 



FIG. 1: Subdivision of a two leg ladder 
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The Hamiltonian of the whole ladder 

H = H^Nxip)) + H 2 {N 2 {p)) + H 12 , (1.7) 

contains in addition an interaction term, which is 
mediated via the two (dashed) links connecting the 
subladders. Each bond in FigQ corresponds to the 
usual t — J Hamiltonian, composing of an electron 
hopping term t(c\ a Cj tCr + cj a Ci^a) and an exchange 
interaction J (Si ■ Sj — nifij/4). The elimination of 
doubly occupied states has been also imposed. The 
indices i,j refer to nearest neighbour sites, a to 
electron spin and the couplings on legs are defined 
as (*', J'). 

The product ansatz: 

fl>{p,N) = ^(pi,JVi(p))V(p2,iV 2 (p)), (1.8) 

describes a state with charge density p and two spa- 
tially separated phases, the first one with charge 
density p\ and N\(p) sites and the second with 
charge density p 2 and N 2 (p) sites. 

If we estimate the expectation value of the Hamil- 
tonian (|1.7|) with the variational ansatz (|1.8fl . we 
get an upper bound 



s{p) < £l(p) for pi< p< p 2 , 

in terms of the right-hand side of 11.611 . 
Note that the interaction term 
1 



(1.9) 



N 



(i,{p)\H 12 \^{p)), 



does not survive in the thermodynamical limit N — > 
oo. Since H1.9[) is a strict upper bound in the ther- 
modynamical limit for any interval 



Pi 



< 



< 



P2, 



we can conclude that e(p) is a convex function of 
p. The product ansatz l|1.8|) with the two separated 
phases represents the true ground state if the upper 
bound l|1.9fl sharply holds. In this case, we derive 
from (|1.6|l a constant chemical potential: 



de e(p 2 )-e(pi) 

/i = — = = no for pi < p < p 2 , 

dp p 2 - pi 

(1.10) 

which corresponds to a discontinuity in the inverse 
function 



p(p) = 



pi for p — ► /i + 
p 2 for p — > po - 0. 



(1.11) 



It is the purpose of this paper, to demonstrate that 
the generic mechanisms, which lead to gaps and phase 



separations, are intimately related with the formation of 
clusters. On ladder systems the formation of clusters is 
prescribed in a natural way by the ladder geometry. 

It is plausible to start with the simplest clusters, de- 
fined by the rungs. In the limit of vanishing hopping 
parameter t' along the legs, the system decouples into 
a product of rung eigenstates. This limit has been 
studied intensively in the literature under various nota- 
tions like "local rung approximation"™ or "bond opera- 
tor theory"it InRefiS, we studied first order corrections 
in the leg hopping parameter t' and compared perturba- 
tive results with exact diagonalizations on a 2 x 8 ladder 
for parameter values 



i=l, a = J = 0.5, 
t' = 0.1,0.2,0.3, a' = J'/t' = 2.7 



(1.12) 



Concerning the quality of the perturbation expansion 
based on the local rung approximation we found 

• good agreement in the regime < p < 1/2, 

• failure in the regime 1/2 < p < 1 . 

We do not think that higher orders in perturbation ex- 
pansion with local rungs will improve the situation in 
the second regime. Instead we are convinced the starting 
point - i.e. the clusters which define the zeroth order 
perturbation theory - has to be changed. 

In this paper we intend to demonstrate on the two 
leg ladder with anisotropic couplings, how the appropri- 
ate clusters which define the zeroth order perturbation 
theory are to be found. We start in Section [D] with an 
analysis of the exact ground state energy per site e(p) 
on a 2 x 8 ladder [cf. Fig. |2{a),(b)] which turns out to 
be almost linear in the charge density p for an appro- 
priate choice of the system parameters (Eqs.(|OJ(a) and 
(JOJ(b)). As explained above the linearity in p indicates 
phase separation into two clusters with charge densities 

Pi=0 P2 = \ for 0<p<i, (1.13) 

Pi = \ P2 = l for \<P<^- (1-14) 

The dependence of the ground state energy per site l|1.6|) 
on the system parameters is very well reproduced by the 
plaquette ground state energies E^(Q) with plaquette 
charges Q — 2,4. 

In Section IIIII we present the perturbation expansion 
based on plaquette clusters. First order corrections, 
which describe the interaction between neighbouring pla- 
quettes, favour the clustering of plaquettes with charges 
Q = 4 - a first indication of phase separation. 

In Section IIVI we investigate under which circum- 
stances plaquette clusters with odd charges Q = 1 and 
Q = 3 are suppressed energetically. 

In the low doping case p > 3/4 the dominance of Q = 
2 and suppression of Q = 3 plaquettes in the ground 
state explains the emergence of tightly bound hole pairs 
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(RefiiS) for an appropriate choice of the rung and leg 
couplings. Finally, a discussion on our results will be 
presented. 

II. APPROXIMATE LINEAR CHARGE 
DENSITY DEPENDENCE OF THE GROUND 
STATE ENERGY 

We start in Figs. H{a), Gib) with Lanczos results of 
the ground state energy per site in the t — J model with 
anisotropic couplings t, a = J/t (throughout this paper 
we have chosen t = 1), t' , a' — J' jt' for the rungs and 
legs respectively on a 2 x 8 ladder: 

(a) a = 0.5, a' = 4.0, t' = 0,0.1, 1.0, 

(b) a = 4.0, a' = 0.5, t' = 0, 0.1, . . . , 1.0 . 

(2.1) 




P P 



FIG. 2: Lanczos results for energies per site of a TV = 2 x 8 
t — J ladder for the parameters l2.ll fa1.CbL For both shown 
cases t' increases from top to bottom as t' = 0, 0.1, . . . 1.0. 

There is a pronounced difference in the p-dependence 
of the ground state energies in Fig. EJa) and Fig. 0{b), 
respectively. In Fig. EJa) we observe a discontinuity in 
the slope p(p) = de/dp at p = 1/2, which generates 
a plateau in the charge density p{p) as function of the 
chemical potential 

p(p) = - P- <p< p+- (2.2) 

The plateau width 

A = p + - p_ (2.3) 

shrinks with increasing values of t'. Note also that the 
chemical potential p(p) — de/dp vanishes for p > 1/2 
for a specific value of t' = t' (a, a'). This feature will be 
discussed in Section Hvl 

Let us next turn to the ground state energy e(p) in the 
domain (|2.1|l fb1 shown in Fig. [2Jb). At t' = 0, there is 
no discontinuity in the slope - i.e. no plateau in p(p) - at 
p = 1/2. The variation with the leg hopping parameter t' 
is much smaller than in case (a). For t' > 0, e(p, t') is only 
approximately linear in the two subintervals < p < 1/2 
and 1/2 < p < 1 with different slopes p- and p + . The 



difference A = p + — p_ increases with t' . Deviations 
from linearity are convex as predicted by (|1.9fl . 

The product ansatz l|1.8fl describes a system with two 
phases: 

For < p < 1/2 with a ground state energy 

£ (p,t') = e(0,t')(l-2p)+e(l/2,t')2p (2.4) 

there is a phase with charge density p± = in the first 
part of the ladder and a second phase with charge density 
p-y = 1/2 in the second part. Phase separation occurs at 
N 1 (p) = (l-2p)N. 

For 1/2 < p < 1 with a ground state energy 

e(p,t') = 26(1/2,0(1 -p)+e(l,*')(2p-l) (2-5) 

the two phases in the first and second part of the ladder 
have charge densities pi — 1/2 and p 2 — 1, respectively. 
Here the phase separation occurs at Ni(p) = 2(1 — p)N. 

For < p < 1/2 the dependence of e(p,t') TTty on 
the parameters t',a,a only enters via £(1/2, t', a, a ). 
Results on a 2 x 8 ladder are shown for this quantity in 
Figs. Ola) and[3{b) with the parameter choices 

(a) a = 0.4, 0.5, 0.6, a = 4.0, 5.0, 6.0, t' = 0, 0.1, . . . , 1.0 , 

(b) a = 4.0, 5.0, 6.0, a' = 0.4, 0.5, 0.6, t' = 0, 0.1, . . . , 1.0 . 

(2.6) 




/ t'/(a-1.4*a') 

FIG. 3: Correction factor /(1/2, t', a, a') of |£7|> for a 2 x 8 
t — J ladder with parameters (I2.6ll fa1.(b1 and p — 1/2 

It turns out that the whole dependence on the param- 
eters t',a,a': 

e(l/2,t',a,a') = 

±EM%t',a,a!)- f{l/2,t',a, a ') (2.7) 

is correctly reproduced by the (2 x 2)-plaquette ground 
state energies E^(Q, t', a, a') with charge Q — 2 up to a 
correction factor /(1/2, t', a, a'). The physical interpre- 
tation of / is given below. 

As is shown in Appendix |X] [cf. HA13MA18[) ]. we 
have different ground states in the regimes I|2.6|l fa1 and 
(|2.6|) fb1 for plaquettes with charge Q = 2. The ground 
state energies follow from the lowest eigenvalues of the 
3x3 matrix l|A9fl . which are computed in a perturba- 
tion expansion with zeroth order contributions (|A15|I and 
(|AL8)l for the regimes (E3(a) and (|2"U|) fbh 



4 



(a) E$(2,t',a,a') = -\ (j' + VJ' 2 + 16) , 

(b) E$(2,t',a,a') = -]- (a+ V« 2 + 16i' 2 



(2.8) 

In Figs. Ofa) and EJb) we have plotted the correction 
factors in l|2.7(l : 

(a) / (a) (l/2,t',a,a') = / (o) (l/2,j' = aV) , 

(b) / (6) (l/2,t',a,a') = / w (l/2,t7(a-1.4a'))- 

(2.9) 

versus the "scaling variables" J' = a'i' and t'/(a— 1.4a'), 
respectively. Note, that all data points for (|2.6|l fa'l and 
(|2.ti|) ('b') almost coincide if we use the scaling variables. 

Let us next discuss the linear p-behaviour in l|2.5[l . 
The dependence on the parameters t' , a, a' enters 
via e(l/2,t',a,a') [^,^(a,),^(b)} and s(l,J = 
a, J' = a't'). Note that for p = 1, e(l, J, J') is just the 
ground state energy per site of the spin 1/2 Heisenberg 
model on a two leg ladder with spin couplings J and J' 
along the rungs and legs, respectively^ The J, J' depen- 
dence 

e(l, J, J') = ^E^(4,J,J')-f(l,J'/J) 

(2.10) 

is correctly reproduced by the plaquette ground state en- 
ergy E^{Q, J, J') with charge Q = 4 [cf. jH| ]: 



Figs. EIa),EIb) and which scales in J'/ J. 



1.25 

S 1.2 
S 

1.15 
, = 11 



• a'=4 
□ a'=5 

♦ ct'=6 



(a) a'=4,5,(5 

tx=0.4,0.5,0.6 




#ffcj a' =0.4,0.5,0.6 
a=4,5,6 



9 12 15 



.17.1 



0.05 0.1 
J'/J 



FIG. 4: Correction factor f(l,t',a,a!) of lETul for a 2 x 8 
i — J ladder with parameters i2.dl (a).(b) and p = 1 



III. PLAQUETTE CLUSTERS IN THE 
GROUND STATE OF THE t - J MODEL ON A 
TWO LEG LADDER 

The results of the numerical analysis in the last section 
motivate us to build up the ground state on the two leg 
ladder from plaquette eigenstates. 



For this purpose the t — J Hamiltonian 

N/4 JV/4-1 

H = t^/ijj(t») + t" h JJ+i( a ") 

3=1 3=1 

(3.1) 

is decomposed into NJ 4 plaquette Hamiltonians [using 
the notation of Ref. I|19|) ] 

hjj(t , a! \a) — 

j [h(4j - 3, 4j - 1, a') + fc(4j - 2, 4j, a')] + 
[h(4j - 3, 4j - 2, a) + h(4j - 1, 4j, a)] (3.2) 

with spin couplings J = tot, J' = t'a' and hopping terms 
t, t' along the rungs and legs, respectively (cf. Fig. [SJ. 

hj,j+x{a") = 
[h{4j - + I, a") + h(ij, 4j + 2, a")} (3.3) 

describes the interaction between neighbouring plaque- 
ttes j,j + 1. This interaction is treated in the follow- 
ing in a perturbative expansion around t" = 0, J" = 0, 
a' = a" = J" It" fixed. In the analysis of the numerical 
results for the ground state energies 1|2.7[) and l|2.1UI) at 
p = 1/2 and p = 1, this interaction generates the correc- 
tion factors f(p = 1/2) and f(p = 1) - shown in Figs. 
|3{a),(b) and^Ja),(b) respectively. 



2 t ,a 4 t ,a 6 



(2.11) 














is shown in 


t,a 




t,a 






Qj=N/4 






t',a' 3 


t",a"~ 


N 



N 



FIG. 5: Decomposition of the two leg t — J ladder into coupled 
(t",a") 4-site plaquette clusters 

The lowest plaquette eigenstates 

h^(t' ia ' ia )^{Q 3 ) = E^(Qj)ip^(Qj) (3.4) 

in the sector with charge Qj = 0,2,4, Qj = 1,3 and 
total plaquette spin and 1/2, respectively, are discussed 
in Appendix^] The eigenstates ip^L (Qj) = ^ {p) (Qj) 
with lowest energy 

EnLoiQj) = E (p) (Qj) y ield the basis 
for the ground state on the two leg ladder in lowest order 
perturbation theory t" = 0. 

To zeroth order in t" (t" = 0) the eigenstates of the 
Hamiltonian (I3. 1|) are given by a product of plaquette 
ground states ip^ p \Qj) 



N/4 

n^ (p) (^) 

3=1 



(3.5) 
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with energies 



N/4 

E n = ^E^iQj). 

i=i 



(3.6) 



The plaquette charges Qj have to add up to the total 
charge: 



Ot 



ot 



(3.7) 



Our analysis in Sec. |n] suggests that only plaquette 
ground states with charges Q = 0, 2, 4: 

V> (p) (0), V (p) (2) for 0<p<l/2 (3.8) 
V> (p) (2), V (p) (4) for l/2<p<l (3.9) 

are involved in the construction of the ground state of the 
two leg ladder. The corresponding ground state energies 
on the ladder: 



E{p,t',a',a,t" = 0) = N^(i)E^'(Q = i) 



i = for 0<p<l/2 
i = 2 for 1/2 < /? < 1 



(3.10) 

+N^(t + 2)E ( ^(Q = i + 2). 

(3.11) 
(3.12) 



are obtained from plaquette ground state energies 
E&'(Q) and the number N^(Q) of plaquettes with 
charge Q 

7V (0) (0) + A r(0) (2) = N {2 \2) + N {2 \4) = N/4 (3.13) 
27V (0) (2) = 2iV (2) (2) +4iV (2) (4) = Q tot (3.14) 

So far we only treated the zeroth order perturbation the- 
ory (|3.1|l {t" — 0). The product states are degenerate 
since the N^(Q) plaquettes with charges Q — 0, Q = 2 
for i = and <3 = 2,Q = 4fori = 2 can be distributed 
in different ways over the ladder (Fig. (5J . 

The interaction energy W(Qj, Qj+i) between neigh- 
bouring plaquettes is derived in Appendix[H]in the frame- 
work of a I s * order perturbation theory in the hopping 
parameter t". The resulting shift AE in the ground state 
energy: 

AE (0) = N {0) (2,2)W{2,2) (3.15) 
for i = 0, < p < 1/2; 

AE {2) = N {2) (2,2)W(2,2) + N ( - 2) (2,4)W(2,4) 

+7V (2) (4,4)1^(4,4) (3.16) 
for i = 2, 1/2 < p < 1. 

can be expressed in terms of the interaction energies 
W(Qj, Qj+i) and numbers N(Qj,Qj + i) of neighbour- 
ing plaquettes with charges Qj,Qj+\. Since W(0,0) = 
W(0,2) = and W(2,2) = -J'/8 JHEJl, Afi(°' is mini- 
mal if AT(°)(2, 2) is maximal: 

(2, 2) = 7V(0) ( 2 ) _ i = E p _ i > (3.i 7 ) 



According to (|B8|) : 

W(2,2) + W(4,4)-2W(2,4) < 0. (3.18) 

AE ( V is minimal if A( 2 >(4,4) (and thereby N^(2,2)) 
are maximal. 

A^ (2) (4,4) = 7V (2) (4) _ i = y(p- 1/2) - 1, (3.19) 
N^(2,2) = N™(2) - 1 = y(l -p) - 1, (3.20) 



7V (2) (2,4) = l 



(3.21) 



The perturbative result of the ground state energy per 
site in the regime Q2.6[l fa) 

s(p) = (l-p)i(V p >(2) + ^(2,2)) 

+\{p- 1/2) (E^ (4) + W(4,4)) (3.22) 

predicts a level crossing in all charge sectors with 1/2 < 
p < 1 if 

£ (p) (4) (2) + Vy(4,4) -W(2,2) = 0. 

(3.23) 

In the regime <|2.6[l fa') the left-hand side with l|2.11|l . 
(123 (a) and 

W(4,4)- W{2,2) = (3.24) 
8 

only depends on J, J'. The solution defines a curve 
J' = J'(J) in the parameter space, where the ground 
state energy (|3.22ll becomes independent of p, the corre- 
sponding chemical potential: 

de 

fj,= — = for 1/2<o<1 and J' = J'(J) 
dp 

(3.25) 

turns out to be zero here. 

Moreover, the ground state energy along the curve 
J' (J) is predicted to be 



s(p,j' = j'(J)) = \ E ^{2)-^ 2 



(3.26) 



Results for J' (J) and e(p. J' — J' (J)) are given in Table 



J = a 


0.4 


0.5 


0.6 


J' = J'(J) 


0.94895 


0.88839 


0.81255 


e(p, J'(J)) 


-0.66215 


-0.65099 


-0.63717 



TABLE I: Results for J' (J) and e(p, J' = J' (J)) for run g spin 
exchange couplings J — a — 0.4, 0.5, 0.6 on the basis of I3.23H 
and EJHl. 
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r+ 




FIG. 6: Lanczos results for the crossing of energy levels for 
P = Q/N > 1/2 and a = 0.5, a' = 4.0(a), a' = 2.7(b) for a 
N = 2x8t — J ladder 



As an illustration, we present in Fig. Elthe ground state 
energy per site e(p,t',a') = E(N,Q,t',a',a = 1/2) /N 
for the charges Q = 8, 10, . . . , 16 at a = 0.5, a' = 4.0 (a) 
and a' = 2.7 (b). 

The crossing of these energy levels at 

, J'(J = 0.5) _ f 0.222 for a' = 4.0 
a' ~ 1 0.329 for a' = 2.7 

(3.27) 

is predicted to change with a' if we keep a = J = 0.5 
fixed. On the other hand, the corresponding ground state 
energy per site 

e{p, J'{J = 0.5)) = -0.65099 (3.28) 

is independent of a' I 

Both predictions (|3.27|) and l|3.28|) are clearly visible 
in the numerical results on a 2 x 8 ladder. 

We also looked for level crossings <|3.23[1 in the regime 



(12.611 (b). It turns out that they occur at t' values t" = 
if > 1, where the perturbative approach is not reliable. 

IV. THE SPECIAL ROLE OF PLAQUETTE 
CLUSTERS WITH CHARGES Q = 2 AND Q = 4 

The variational ansatz (|3.5I) with plaquette clusters 
only involves cluster eigenstates ()3.8|l . Ij3.9j) with even 
charges Q — 0,2,4. Such an ansatz makes sense if pla- 
quette clusters with odd charges are suppressed energet- 
ically: 

Ai = £^(0) + £ (p) (2)-2£^(l) <0, (4.1) 

A 3 = £ (rt (2) + £ (p) (4)-2£ (p) (3)<0. (4.2) 

In Appendix ^ we discuss the low-lying eigenstates of 
plaquette clusters with charges Q = 0,1,2,3,4. The 
ground state energies for Q = 0,1,4 are unique in the 
sense that there are no level crossings by variation of 
the parameters t' , a' , a. The ground state energy of the 
Q = 2 cluster is given by (|2.8|l (a) and 12.8|l (b). respec- 
tively. 

It turns out that the inequalities (|4.1|) and (|4.2|l are 
indeed satisfied in the regimes l|2.1(l (a) and l|2. 111 (b). re- 
spectively. 

The suppression of Q — 3 plaquettes in the regime Ij4.2|l 
has immediate consequences for the mobility and correla- 
tions of holes as they are discussed by Siller, Troyer, Rice 
and Whitoi& for the low doping case (p > 3/4). In the 
regime A 3 < hole pairs are confined in Q — 2 plaquettes 
and cannot move in the antiferromagnetic background of 
the Q — 4 plaquettes. This can be seen in the first order 
perturbation theory (cf. Appendix [B| for the interac- 
tion between neighbouring plaquettes, which forbids the 
transition 

Q] = 2 Q j+ i = 4 - Q;. = 4 Q' j+1 = 2 (4.3) 

with charge transfer AQ = 2. Therefore the hopping of 
the hole pair confined in the Q — 2 plaquette is sup- 
pressed in the antiferromagnetic background. 

On the other hand hopping of a single hole contained 
in a plaquette with Q = 3 is possible since the transition 

Qj = 3 Q J+ i = 4 - Q'j = 4 Q' j+1 = 3 (4.4) 

is not forbidden. 

It is interesting to study the energy difference A3 (a, a') 
for the a, a' values of Ref. (fish . There the leg coupling 
a' = J' /t' has been chosen to a' — 0.35 whereas the rung 
coupling varies a > a' = 0.35. The hopping parameters 
are equal t — t' . For large a > 4 this corresponds to our 
regime (I2.6|l (b). 

In Fig. [7| we have plotted the difference A3 (a, a' = 
0.35). For comparison we have also included the differ- 
ence A 3 (a = a') with symmetric couplings along the 
rungs and legs, respectively. The latter drops monotoni- 
cally with a — a' and has zero at 

a = a' = = 0.436... (4.5) 

v 21 
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FIG. 7: Energy differences A 3 (a,a') for (1) a' = 0.35 and 
(2) a = a' (in both cases: t = t' = 1). 

In contrast the difference A3 (a, a' = 0.35) in the asym- 
metric case a > a' = 0.35 first increases and has a flat 
maximum at a ~ 0.43 and then drops with a zero at 

A 3 (a ~ 0.985, a' = 0.35) = 0. (4.6) 

The zeros ()4.5|l and (|4.6(l of A3 mark the transition where 
Q = 3 plaquettes are substituted by Q — 2 plaquettes 
and 2 holes combine to a pair. We expect that the hole- 
hole correlation length has a maximum at this transition 
point. Indeed there is a maximum of the hole- hole corre- 
lation length at a = 1.2, as determined in Ref. I|18r) from 
a DMRG calculation on a 40 x 2 ladder. It might be an 
accident, that the value a = 1.2 is quoted as well as a 
lower bound for phase separation in the 2D t — J model 
with isotropic couplings (Putikka et al. 5 ). We therefore 
calculated (|4.2() as well for the isotropic case a = a', 
where the zero 1(4. 5[) is formed quite below a = 1.2. 

It would be interesting to see whether the maximum of 
the correlation length is shifted as well to a smaller value 
in the symmetric case a — a' . 

We are aware of the fact, that our considerations in the 
low doping regime (p > 3/4) are based on the product 
ansatz (|3.5|l with plaquette clusters with charges 

Q = 2,4 for A 3 < (4.7) 
Q = 3,4 for A 3 > (4.8) 

The interaction between neighbouring plaquettes is ne- 
glected. This is justified for (|4.7I) . if A3 is sufficiently 
negative, as was demonstrated in Sees. ITT1 and ITTT1 On 
the other hand, these interactions cannot be neglected 
in the regime l|4.8fl . where first order perturbation the- 
ory allows the hopping (|4.4|l of Q = 3 plaquettes in the 
antiferromagnetic background of Q — 4 plaquettes. In 
this case, first order perturbation theory leads to an ef- 
fective Hamiltonian on a chain with nearest neighbour 
couplings. The effective degrees of freedom at each site 
and their nearest neighbour interactions are defined by 
the ground states of the Q = 3 and Q = 4 plaquettes. 
We expect, that such an effective Hamiltonian will induce 
a convex curvature in the p-dependence of the ground 



state energy, which would indicate that we are beyond 
phase separated phase. We observed this curvature (for 
< p < 1/2) with an effective Hamiltonian based on 
rung clusters and their interactions^ 



V. DISCUSSION AND PERSPECTIVES 

The formation of clusters in the ground state of a quasi 
one-dimensional system has important consequences for 
its physical properties, e.g. the phase diagram at zero 
temperature. 

In case of the t— J model on a two leg ladder with asym- 
metric couplings (a — J/t,a' = J'/t'; regime H2.6|) (a~) 
a /a' <C 1 and regime ((2.6(1 fb) a /a' ^S> 1) plaquette clus- 
ters with even charges Q = 0, 2, 4 play the dominant role 
and explain the charge density and (a, oe',t') dependence 
of the ground state energies, as was demonstrated in Sees. 
Illll and llVl First order perturbation theory for the inter- 
action of neighbouring plaquettes favour the condensa- 
tion of plaquettes with the same charge (Q = 2,Q = 4), 
which can be interpreted as a signal for phase separation. 
Of course, this can happen only if plaquette clusters with 
odd charges Q = 1, Q = 3 are suppressed energetically, 
which means that the energy combinations i|4.1|) 14.2(1 are 
sufficiently negative. 

The low doping regime (p > 3/4) is of special interest. 
It has been demonstrated in Ref. flSr) that the hole-hole 
correlations are small in the t — J model with asymmetric 
couplings on legs (a' = 0.35) and rungs (a > 4). In this 
regime the ground state is very well described with pla- 
quette clusters of charge Q = 2, Q = 4. The Q — 4 pla- 
quettes generate the antiferromagnetic background. In 
each Q = 2 plaquette a pair of holes is confined. The 
holes can be deconfined only if the resulting two plaque- 
ttes with charge Q = 3 are energetically preferred [cf. 
(|4.2|) for A 3 > 0] . As was demonstrated in Sec. |TV| this 
happens for smaller a values (a' — 0.35 < a < 0.98 in the 
asymmetric, a = a' < 2/y/21 in the symmetric case). In 
this regime the ground state is more complex [Ref. ijlSI) ]. 

In our approach based on a product ansatz with clus- 
ters we have two possibilities to exploit the ground state 
in this regime: (A) We compute the effective Hamilto- 
nian which describes the interaction between Q = 3 and 
Q = 4 plaquettes perturbatively. (B) We improve the 
quality of the product ansatz with larger clusters such 
that hole-hole correlations at larger distances are prop- 
erly taken into account as well. 
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APPENDIX A: EIGENSTATES OF THE 4-SITE 
CLUSTER 



a) Q = 0, N c (0) = 1 



The Hamiltonian of the 4-site cluster (cf. Fig. |HJ is 
defined in i|3.1 (I (for j = 1). 



2_ f,a' 4 



t,a 



Qi 



t,a 



(1,2,3,4) 



i t',a' 3 



FIG. 8: Notation and shape of the used 4-site cluster - build- 
ing block for the two leg t — J ladder 

We are now going to construct all N C (Q) eigenstates 
for cluster charges <2 = Qi=0, 1,2,3,4- assuming equal 
numbers of spin-up (+) and spin-down (— ) charges for 
even Q and an excess of one spin-up particle for odd 
cluster charges Q. Moreover, we will assume t,J,t',J' > 
0. 

The difference in coupling parameters for rungs (t, J) 
and legs (t',J') subdivides the basis states for each Q 
into the classes {\III)} and {\I),\II)} 



\III) 



Q = l 


Q 


= 2 


Q = 3 


\ 







+ 


+ »/• 


+ 






+ o\ 


+ 


-)■ 




0/' 







+ 



(Al) 



and^ 



\II): 



- 
+ 

+ 

- 











(A2) 



(A3) 

concerning the behaviour with respect to transformations 
R: (1,2,3,4) — > (1,3,2,4), (A4) 
i.e. the interchange of legs and rungs. 



£(P)(Q = 0) = 0. 



(A5) 



b) Q = l, N c (l) = 4 

We introduce \n) as Q = 1-basis state with the 
spin-up particle (+) at plaquette position n (\n) = 
c+ T |0». 

Using the basis: 

|1,t,C) = i(|l)+r|4)) + C(|2)+r|3» (A6) 
with t = ±1,C = ±1, we obtain the eigenvalues: 
EW(Q=1,T,0 = -C(t + rf). (A7) 

c) Q = 2, N c (2) = 12 

Introducing |m, n) as Q = 2-basis state with spin- 
up particles (+) at site m and (— ) at position n 
(\m,n) = c+jC^JO)), we first compute the action 
of the t — J Hamiltonian onto the singlet (r = 1) 
and triplet states (r = —1): 

|l,r,C) = |l,2) + r|2,l)+C(|3,4)+r|4,3» 
|2,r,C) = |2,4)+r|4,2)+C(|l,3)+r|3,l)) (A8) 
|3,r,C) - |l,4)+r|4,l)+C(|3,2)+r|2,3)) . 

In this basis the action of the t — J Hamiltonian 
results in a 3 x 3 matrix for the singlet (r = 1) 
sector 



/ -J 
—J' 

V-(i+o*' -(i+o* 




(A9) 



The ground state energy E in the singlet sector 
(t = 1 with C = 1) is found from the solution of 
the third order equation 

- (J + E){J' + E)E + i(E + J)t 2 + A{E + J')t' 2 

= (A10) 

which can be easily solved in the symmetric case 
t' = t, J' = J: 



E 



(j+ VJ 2 + 32) . (All) 



In the asymmetric case H2.6[l fal with a — J/t <C 
a' = J' jt' one can derive an iterative solution treat- 
ing the term 



{J' + E)E + At 2 = -4 



E + J' 
E + J 



(A12) 
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on the right-hand side of l|A12(l as a perturbation. 
The resulting ground state energy in first order of 
this perturbation reads: 



AE n = 



-4 



e (p) 



AE a 
J' 



E (v) 



J V J' 2 + 16t 2 



where 



E (p) 



J' + v J' 2 + 16t- 



(A13) 
(A14) 

(A15) 



E.g. for J = 0.5, J' = 4, t' = 1 the correction 
term (|A14|) yields a 3% contribution to the ground 

state energy E, such that the zeroth order Ea is 
already a very good approximation. 

In the regime (|!H)|)(b ) with a' = J'/t' < a = J/t 
we derive from (jAlOfl a corresponding ground state 
energy in a first order perturbation: 



E, 



(p) 



E. 



(p) 



AE h 



AE h 



A E b + J 



E b + J' V J 2 + 16£ /2 



where 



(A16) 
(A17) 

(A18) 



Let us now turn to the triplet sectors (t = —1). 
Here the 6 eigenstates and corresponding eigenval- 
ues turn out to be 

|1,(-1,-1)) E = 

|2,(-1,-1)) E = 

|1,(-1,1)>±|3,(-1,1)) E = ±2t' 

|2,(-1,1))±|3,(-1,-1)) E = ±2t. 

d) Q = 3, AT C (3) = 12 

We introduce |m, n) as Q = 3-basis state with the 
hole (0) at plaquette position m and spin-down 
electron (— ) at position n. The creation opera- 
tors of the two spin-up electrons (+) and that of 
the spin-down electron (— ) are ordered according 
to the increasing site number (cf. Fig. |SJ)A. E.g. 
1 1,3) = c + c +c+|0). 

In the basis: 



|1,t,C) = (|4,2> + t|1,3»+C(|3,1)+t|2,4» 
|2,t,C) = (|4,3)+r|l,2))+C(|3,4)+r|2,l)) 
|3,t,C) = (|4,1> + t|1,4»+C(|3,2)+t|2,3» 

(A19) 

the t — J Hamiltonian reduces to the following 3x3 
matrix for the eigenvalues with x = J/2, x' = J'/2: 



/ Crf -x -C* • 
-£t-x' Crtf- 
y -(t + x C,Tt' +x' -x- 



(A20) 



The eigenvalues are given by: 



£(p)(t,C) = -(x + x')± V^W+CSW, 

C(rt' - t) (A21) 

with 

A (T » = (x 2 -xx 1 + x' 2 ) + {t 2 +Ttt' + t' 2 ), 
B {t) = -{2t + Tt')x + {2rt' + t)x' . 

(A22) 

The 4 eigenstates corresponding to the eigenvalues 
C(Tt' — t) which are independent of the spin cou- 
plings J, J' have maximal total spin S = 3/2. The 
remaining ones have total spin S = 1/2. 

For the choices (|2.1f) fa.b) of plaquette parameters 
(t, J, t', J') the ground state is given by one of the 
J, J'-dependent (r = l)-states with energy: 



E(f)(r = l,C) = -(x + x' + v^ (1) + (B(K 



(A23) 

The ground state in the symmetric case t' = t, 
J' = J 




is twofold degenerate with respect to the quantum 
number £ = ±1. This degeneracy is lifted in the 
asymmetric case l|A23|l if the term B^ 1 ' l|A22(l is 
nonvanishing. 

e) Q = 4, A^ c (4) = 6 

We introduce |m, n) as Q = 4-basis state with the 
two spin-up particles (+) at plaquette positions m 
and n. Again, we use the definition that all creation 
operators of the four electrons 2(+),2(— ) act in 
order of increasing site number (cf. Fig. |SJ on the 
vacuum |0) (e.g. |1,3) 

In the basis: 



c i T C 2~l C 3~T c 4| 1 0) ) ■ 



|l,r) - |l,3)+r|2,4) 

\2,t) = |l,2)+r|3,4) (A25) 

|3,r) = |l,4)+r|2,3) 

the t — J Hamiltonian reduces to the matrix 



-J 



AtJ 

-J' AtJ' 
\AtJ AtJ' -(J+J') 



(A26) 



with At = (1 + r)/2 = 1,0. 

The eigenvalues for r = — 1 arc simply given by: 



-J, - J', -{J+J'), (A27) 
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whereas the case r = 1 yields: 
£(f)(Q = 4,+) = 



-(J + J') ± \J J 2 + J' 2 — JJ', . (A28) 

The ground state for all parameter values is 
uniquely given in the r = 1 sector. 



is forbidden in first order perturbation theory. 
Therefore the nonvanishing matrix elements: 



--W{Q h Q j+l )5 Q , j>Qj 5 Q , j+i>Qj+1 (B5) 



APPENDIX B: FIRST ORDER PERTURBATION 
THEORY 



are necessarily diagonal for the pairs of interest 



To see the effects of first order perturbation theory in 
the hopping parameter t", we start from the transition 
matrix elements: 



N/4 
3 = 1 



N/4 \ 



IN/4 

i"(ri4 p) (^ 

\i=i 

N/4 

= t"Y l A j<j+1 (4"\Q' j )4"\Q' j+1 ) 



(«') 



with 



(Bl) 



(B2) 



The interaction Hamiltonian hj j+i (a') between neigh- 
bouring plaquettes is illustrated in Fig. |5J The hopping 

part hjj +1 is active if the two links {xx') (yy') are oc- 
cupied by one electron and one hole, respectively. This 
means in terms of occupation numbers nix), n{x') on the 
sites x, x' n(x) = l,n(x') = or n(x) = 0,n(x') = 1. 
Therefore, the hopping term induces a charge exchange 
by one unit 



■,(*) . 



(Qj + l,Q j+1 -l) (B3) 



between neighbouring plaquettes. Note in particular, 
that charge exchange by two and more units 

(0, 2) «-> (2, 0) (2, 4) (4, 2)(Q 3 - + 1, - 1) 

(B4) 



(Qi,Qj+i) = (0,2), (0,4), (2, 4), (0,0), (2, 2), (4, 4). 

(B6) 





FIG. 9: Leg couplings (a;, a;'), (y,y') linking the neighbouring 
plaquettes j and j + 1 on the ladder. 



They arise from the spin exchange part h>j)+i which is 
active if both sites x, x' (cf. Fig. are occupied with 
one electron: 



• (J) 

L 3,j+1 



2S{x)S{x') + - ) n(x)n(x') . (B7) 



Here S(x) and S(x') are spin operators at sites x and x' 
and we get for the interaction energies W(Qj, Qj+i) for 
plaquette pairs (|B6|I : 



U = 



,// /QjQj+i 

—t a 

32 



(B8) 



if Qj and Qj+i come from Ea. (|B6 
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